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We model spin blockade for optically excited electrons and holes in a charged semiconductor 
quantum dot. We study the case where the quantum dot is initially charged with a single electron 
and is then filled with an additional, optically excited electron-hole pair, thus forming a charged 
exciton (trion). To make contact with recent experiments, we model an optical pump-probe setup, 
in which the two lowest quantum dot levels (s and p shells) are photoexcited. Using the Lindblad 
master equation, we calculate the differential transmission spectrum as a function of the pump-probe 
time delay. Taking into account both spin conserving and spin- flip intraband relaxation processes, 
we find that the presence of the ground-state electron spin leads to an optical spin blockade at 
short delay times which is visible as a crossover between two exponential decays of the differential 
transmission. To make predictions for future experiments, we also study the dependence of the 
spin-blockade on an external magnetic field. 



I. INTRODUCTION 

One of the promising solid-state implementations for 
the realization of quantum computing that has been un- 
der intense study over the past years involves the use 
of a single electron spin confined to a charged quantum 
dot (QD)^ The discrete QD enersv structure allows for 
long spin lifetimes, e.g. exceeding one second in electri- 
cally defined GaAs QDs^, in comparison with the bulk 
materials or semiconductor nanostructures of higher di- 
mension. However, it has been shown that the inhomoge- 
neous dephasing time T2* in GaAs QDs in the presence of 
an unpolarized ensemble of nuclear spins in the QD ma- 
terial is of the order of ~ 10 ns, while the intrinsic spin 
coherence time T2 can reach values beyond 1 ji^. The de- 
coherence time is relevant for quantum information appli- 
cations where it should exceed the elementary quantum 
gate operation time by a substantial factor. Electrical 
control of single spins has been realized in timescales of 
about 50 to 100 nanosecond while ultrafast optical 
pulses have been shown to allow ensemble spin manipu- 
lation in picosecond timescales^ ^ and arbitrary coherent 
single-spin rotationd^. 

Spin blockade, more generally known as Pauli block- 
ade, describes a situation where an electronic process is 
inhibited for certain spin configurations because the ener- 
getically accessible final states are forbidden by the Pauli 
exclusion principle (Fig. [T]). In the electric transport be- 
tween coupled quantum dots, spin blockade can prevent 
an electron to access an energetically favorable path due 
to spin conservatiorP (Fig.jl^). Observations of the leak- 
age current in the spin-blocking regime have allowed the 
study of spin decohere nce m echanisms and in particular 
the role of nuclear spin^^^'^. 

In optical experiments, Pauli blocking effects, also 
known as phase space filling^^, are commonly observed 
in absorption spectra at high photoexcitation intensities 
or when ground state carriers are present. Spin block- 
ade of the lower Zeeman branch in a singly charged QD 
in strong magnetic fields has been studiecP^, and in a re- 
cent pump-probe experiment'"^^, signatures of optical spin 




FIG. 1. (a) Spin blockade in charge transport through a 
double QD, connected to source (S) and drain (D) leads. A 
parallel spin configuration (spin triplet) can lead to a blocking 
of the current, (b) The optical analogy of spin blockade in 
the intraband relaxation between two QD levels s and p in a 
pump-probe set-up. 



blockade have been observed in the transmission spec- 
tra of a charged QD. Lifting of spin blockade is typi- 
cally more difficult to observe for optically excited car- 
riers than for transport setups because electron-hole re- 
combination processes can be much faster than the spin 
coherence and relaxation times. 

In this paper, we describe an optical spin blockade ef- 
fect in a charged quantum dot with two photoexcited 
energy levels that play the role of the two coupled QDs 
in transport. (Fig. ^p) We show the signature of in- 
terlevel (intraband) spin relaxation on the differential 
transmission signal in a pump-probe setup (Fig. |2| and 
draw the analogies between optical and transport exper- 
iments. For small QDs, the admixture mechanisms due 
to spin-orbit coupling play a smaller role^^ and direct 
spin-phonon mechanisms need to be taken into account. 



II. THEORETICAL MODEL 

We study a quantum dot in a cubic semiconductor (e.g. 
GaAs) charged with a single electron. For self- assembled 
quantum dots, lateral dimensions are significantly larger 
than their height, and we thus we assume a circular quan- 
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turn dot in a parabolic confinement potential character- 
ized by a frequency 00. In analogy with atoms, single par- 
ticle eigenstates in QDs are typically labeled as s, p, d,... 
shell, which for our model correspond to n = 0,1,2,... 
harmonic oscillator states, with n = Ux -\- riy the total 
quantum number. Including spin, single particle states 
in the conduction band are degenerate with respect to 
spin Jz = ±1/2 in the absence of a magnetic field (for 
circular QDs). In the valence band, heavy hole (total an- 
gular momentum Jz = ±3/2) and light hole {Jz = ±1/2) 
states are split due to confinement by an energy Ai^. 
Here, we will consider only heavy hole states, assuming 
that the heavy hole-light hole mixing near the band edges 
can be neglected. 

The system is optically excited by a strong pump pulse 
that is resonant to the first excited QD state (p-shell), 
creating an electron-hole (e-h) pair with specific angu- 
lar momentum depending on the pulse polarization. Ac- 
cording to the optical selection rules, a right- (left-) cir- 
cularly polarized a± pulse excites a Jz = Tl/2 elec- 
tron and Si Jz = ±3/2 hole, creating an excited trion 
state (Fig. [2|. Depending on the spin polarization of the 
electrons, the (sp) trion state can be an electron singlet 
(total trion angular momentum Jz = ±3/2) or triplet 
{Jz = ±5/2,±3/2,±l/2)i^ The singlet and triplet 
states are split by an energy Age due to electron-electron 
exchange interactions, which is typically of the order of a 
few meV. In our model, we assume that the pump pulse 
width is much broader than the singlet-triplet splitting 
Age and thus the latter can be ignored. 

We focus on th interlevel relaxation of the photoex- 
cited electron, i.e. relaxation from the (sp) trion state to 
the (ss) trion (see Fig. [2|. Since the latter can only be 
an electron singlet, the relaxation rate depends strongly 
on the excited trion state. If it is a spin singlet, inter- 
level relaxation takes place through phonon emission on 
a timescale of a few tens of ns. On the other hand, if 
it is a spin triplet, a spin flip mechanism is required for 
the relaxation to take place. This will typically involve 
spin-orbit coupling in combination with phonon emission 
and take a much longer time as compared to the spin- 
conserving relaxation. 

In our model, we use the HamiltoniarP^ 



H = Hq -\- Hl ± He , 



where 



ncr na 

describes non-interacting electrons and holes and 
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is the coupling to the optical field, where ej^^ (^na)' ^na 
{hna) are the creation and annihilation operators of an 
electron (hole) in the n-th quantum dot level (n = s^p) 
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FIG. 2. Energy levels of a charged quantum dot under res- 
onant photoexcitation with right- (left-) circularly (cr±) po- 
larized light of the first excited level (p shell), indicated by 
curved blue arrows. The electron (hole) spin in the lowest QD 
level (s shell) is denoted by -l (TT, while t*, denotes an 
electron spin in the excited QD level (p shell) . Straight green 
arrows indicate relaxation processes with Fc the intraband 
spin-conserving relaxation rate, Fg the intraband spin-flipping 
rate, and Frad the interband radiative recombination rate. 



with spin cr = ±^(cr = ±|), E^^ (^na) siugle- 
particle energies of the QD levels, d the interband dipole 
moment, and E{t) the electric field. For a more compact 
notation, we use the notation a =t, | when a =|, t- 

The last term in the Hamiltonian Eq. ([T]) describes 
Coulomb interactions. 
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where only terms that conserve the number of particles in 
each QD level are included. This is a reasonable approx- 
imation for very small QDs in which interlevel spacing is 
much larger than Coulomb interaction. Such terms lead 
to density dependent energy shifts, as we will discuss in 
the next section. 

Intraband relaxation of electrons from the p to the s 
shell is described with the Lindblad operators 



st 



(5) 



Similarly, hole relaxation is described by the operator 



h 



(6) 



The dynamics of the density matrix p describing the elec- 
tronic state of the quantum dot is given by a master 
equation in the Lindblad form (h = 1 throughout the 
paper). 
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where r = e,h and 



To 
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are phenomenological electron intraband spin-conserving 
and spin-flip relaxation rates. In a recent pump-probe 
experiment on a CdSe/ZnSe quantum dot^^, the two re- 
laxation rates have been estimated to be of the order 
of Fg ~ 0.01 ps~^ and Fc ~ 0.1 ps~^ respectively, corre- 
sponding to two well separated time scales. 

Hole spin relaxation has been found to be much slower, 
of the order of = 1/F^^ - 20 ns^^ in both CdSe and 
In As quantum dots and can be safely ignored here. We 
will only consider hole charge relaxation Fh = F^^. 



III. EQUATIONS OF MOTION 

To compare with pump-probe experiments, we calcu- 
late the differential transmission signal AT/T 
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where Ton (Toff) is the probe pulse transmission coeffi- 
cient when the pump pulse is on (off), and P(3) is the 
induced polarization in frequency space in third order in 
the optical fleld. 

The polarization is connected with the off-diagonal 
density matrix elements 



dY^ P 



(10) 



where 



(Pna) = {hnaena) = TT[hnaenap] (H) 



describes the interband excitation of an e-h pair in level n 
with spins a and a respectively. Here, we have introduced 
the average (• • •) = Tr[---p]. Using Eq. ^ with only 
two QD levels (s and p) per band and factorizing all 
four-operator expectation values within the Hartree-Fock 
approximation, the polarization dynamics is described by 



-dE{t){l - N^l, - Nl) 
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where 



(13) 



are electron and hole populations, and we have defined 
t^nm = = ^n'^^-Kt- The primed summation 

runs over all states {ma'} ^ {^cr}, and polarization de- 
phasing is described with a phenomenological dephasing 
rate 7p. 

The flrst three terms of Eq. (fl2| correspond to the 
semiconductor Bloch equations^^. The second term is 
the standard phase space fllling term due to Pauli block- 
ing, while the third term describes the renormalization of 
single particle energies due to Coulomb interactions. The 
last two terms describe a population-dependent dephas- 
ing of polarization due to electron and hole relaxation. 

The dynamics of electron and hole populations is de- 
scribed by similar equations of motion. 



(Tia2 

+ iSrhThNl^A^ - N^,){5„s - Snp), (14) 

with r = e^h and the population relaxation rate. 
Again, the last two lines in Eq. (14) describe the effect 



of intraband p ^ s shell relaxation. 

Since in pump-probe experiments the measurable 
quantities are at least third order in the optical fleld, 
the above equations may be expanded in terms of in- 
creasing order in E{t), i.e. P^a = PnJ + PnJ + 0{E^) 
and 7V;, = Nn^^^ + K^^^ + 0{E^). Note that K^^^ 
is essentially the ground state population, which van- 
ishes for undoped systems. In our case, assuming that 
the ground state electron lies in the lowest QD level, 

iVna^^ = = ^re^ns^ta whcrC is the S-shcll fllliug 

factor. 

In this manner we obtain a closed set of equations up 
to third order in the optical fleld, which are written ex- 
plicitly in Appendix A. In the next section we will discuss 
their analytical and numerical solutions and calculate the 
differential transmission signal. 



IV. RESULTS AND DISCUSSION 

A. Analytical Solutions 

The equations derived in the previous section can now 
be solved numerically for any exciting laser fleld E{t). 
In the special case of ultrashort pump and probe pulses 
that can be described by delta functions, Eqs.( Al )-( A5 ) 
can be solved analytically. Even though in this case all 
QD levels can be excited (which is not the case in the 
experiment), analytical expressions provide useful insight 
for the dynamics and we will discuss them briefly in this 
section. 

We assume an optical fleld that consists of two laser 
pulses propagating with time delay r with respect to each 
other, i.e. it has the following form (at the QD): 



E{t) = T^probe(^) + ^pump(^ + T ) 



(15) 
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FIG. 3. Schematic representation of a typical pump-probe 
setup. The system is photoexcited by a strong pump pulse 
followed by a weaker probe pulse after time delay r. The 
signal emitted in the direction of the probe pulse is measured 
as a function of r. 



where Ei{t) = El5{t), i = pump, probe, and E^^^^^ 
j^^pump^ is the amplitude of the probe (pump) pulse that 
arrives at the system at time t = {t = — r) (Fig. |3|. 

Using Eq. (15) in the equations of motion (for details 
see Appendix A), we obtain the interband polarization 
in first order in the optical field. 



P£)(t) = zd{l - [Er''''e-^^-^'e-^-^'0{t) 

+^P^"^Pe-'^-('+")e-^-('+")^(t + r)j (16) 

which consists of two parts due to the two pulses in the 
optical field. 

For quantities that are second or third order in the 
optical field, we will only retain terms that are up to 
first order in the probe pulse, assuming that it is much 

In^this case, 
has the 



weaker than the pump (£^q^°^® <C Eq^^ 
the solution for the hole population, Eq. iAA) 
form 



Nt{t) = \d\\l-u^^^)Er-^ 

X [0{T)e-^^''0{t) + 0{-T)e-^^'^'+^^0{t + r)] } (17) 

which describes the creation of hole population in the n- 
th shell either from the pump pulse only, or from both 
the pump and probe pulses. Here we defined = 

For the electronic populations we obtain similar ex- 
pressions, but 7Ar is replaced by a level-dependent relax- 
ation rate = 7a/" + ^npiKa^c + Ka^s) and there are 
additional terms of the form 



that describe the rise of the 5-shell electron population 
due to interlevel relaxation. These terms also appear in 
the solution for the third order terms PAo- , and lead for 
a spin-dependent increase of the differential transmission 
signal as a function of the time delay. The exact expres- 
are included in Appendix B. 




FIG. 4. Differential transmission signal AT/T (in arbitrary 
units) as a function of time delay r and probe pulse energy huu 



for an unpolarized ground state electron, z/J^ 



0.5, and 



sions for Nn^"^ and Pna 



linearly polarized, Gaussian pulses with duration Tpump = 700 
fs, Tprobe = 180 fs. For this plot, we have used the parameters 
Es 2110 meV, Ep 2210 meV, 7p = 5 ps, Fh 0.1 ps, 
Fc — 15 ps, Fs — 170 ps, 7Ar — 480 ps, and Unm — 0. 



B. Zero magnetic field 

In this section, we discuss the results from our nu- 
merical calculations of the differential transmission signal 
for Gaussian pulses similar to the experiment of Ref. iMl 
Fig.|4]shows the imaginary part of the nonlinear polariza- 
tion P^^\uj) for the case of an unpolarized ground state 
electron as a function of the time delay r between the 
pump and probe pulse and the probe pulse energy. There 
is a single peak at the 5-shell trion energy Es that in- 
creases with time delay for tens of picoseconds, in agree- 
ment with the experimental findings. This slow increase 
of the signal is a signature of intraband relaxation from 
the p to the s shell, for which spin-conserving and spin- 
flipping mechanisms contribute, since the optical pulses 
are linearly polarized and the ground state electron un- 
polarized. 

A more detailed description of the dynamics is shown 
in Fig. [5] which depicts snapshots of the signal for specific 
time delays. For r = — 2 ps, for which the probe pulse 
precedes the pump, there is a small signal that arises 
from the interference between the two pulses and is char- 
acterized by oscillations with frequency Ep — Eg. 

For r = 2 ps, when the probe pulse arrives right after 
the pump, the situation is different as the pump pulse 
has created an e-h pair in the p-shell. The hole relaxes 
almost immediately to the s shell in the valence band, 
and as the probe pulse arrives, it can either recombine 
with the ground state electron, or block the probe pulse 
absorption (bleaching) , thus leading to an increase in the 
transmission. 

For r = 20 ps, the electron has relaxed to the 5-shell 
only if it is in the singlet configuration. Given that only 
two of the triplet states are bright, this leads to an addi- 
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FIG. 5. Differential transmission signal as a function of probe 
pulse energy for different time delays r. All parameters are 
as in Fig. [4] The dashed line marks the signal right after 
excitation by the pump pulse, that corresponds to bleaching 
due to hole interlevel relaxation. The dotted line marks the 
expected signal for full electronic relaxation. 
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FIG. 6. Same as Fig. [5] but including Coulomb interactions 
{Unm = 0.1 meV). 



tional increase of the signal by a factor of ~ 1.3 in com- 
parison to the signal at r = 2 ps, that remains constant 
for tens of picoseconds until spin flips can take place. 

In the above results we have ignored the role of 
Coulomb interactions. Their contribution is shown in 
Fig. [6j which shows the differential transmission signal 
for the same parameters as in Fig. [Sj but with the ad- 
ditional terms arising from Coulomb interactions. It is 
clear that their main effect is a shift of the fundamental 
trion resonance for very short timescales, but their role 
is diminished for larger time delays. This is in agreement 
with the results of Ref. |T8l where the role of Coulomb 
correlations has been studied. 

Fig. [Tj^a) shows the effect of electron spin relaxation 
on the differential transmission AT/T for a cr+ pump 
pulse and a a =1 ground state electron (in which case 
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FIG. 7. Differential transmission signal as a function of time 
delay r, (a) with and without the electron spin relaxation and 
(b) with and without hole relaxation for a a =| ground state 
electron and cr+ polarized pump pulse. All other parameters 
Fig.U 
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spin relaxation is necessary for interlevel relaxation) . For 
small delay times r, the slow spin- flip processes do not 
contribute and the signal exhibits a sharp increase due 
to hole relaxation. At larger time scales, the role of spin 
relaxation becomes evident by the slowly increase of the 
signal, the absolute maximum of which depends on the 
spin relaxation rate. As shown by the dotted line in 
Fig.[7|^a), there is no increase of the signal in the absence 
of spin relaxation. 

On the other hand, hole relaxation plays an important 
role at early timescales. This is shown in Fig.[7|^b) where 
the differential transmission signal is plotted for different 
values of T^. In the absense of hole relaxation, the signal 
increases slowly due to electron spin relaxation. 



C. Finite magnetic field 

In the presence of an external magnetic field, more 
spin relaxing mechanisms are allowed, thus enhancing the 
spin- flipping relaxation rate. It has been shown in Ref. [151 
that spin relaxation in QDs is produced by a variety of 
mechanisms that can be separated in two groups: direct 
spin-phonon coupling, and admixture mechanisms due 
to spin-orbit coupling. In both cases though, the finite 
magnetic field leads leads to a ^ 5^ dependence of the 
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FIG. 8. Differential transmission AT/T as a function of 
pump-probe delay time r and magnetic field B for a circu- 
larzly polarized (a) cr+ and (b) cr_ pump pulse. We assume 
the temperature is low enough for the ground state electron 
to be polarized. All other parameters are as in Fig. |4] In 
(a), the large spin-blockade signal for low magnetic fields is 
strongly suppressed at higher fields. 



spin relaxation rate between different orbit als. For the 
quantum dots considered her^^^^ and magnetic fields 
up to 5 T, Zeeman splitting is much smaller (~ fieV)^ 
than the interlevel spacing (50 — 100 meV) and its role 
is insignificant. Thus, the admixture of different spin 
states plays a lesser role and the dominant spin flipping 
mechanism is the direct spin-phonon coupling. 

In Fig. |8| the differential transmission signal at the 5- 
shell resonance is shown as a function of time delay r 
and magnetic field B for right and left circularly polar- 
ized pump pulses. Assuming that the temperature is low 
enough for the ground state electron to be fully polarized 
by the applied magnetic field, a cr_ pulse leads to a well 
defined spin blockade regime, as shown in Fig. jsjb). For 
low magnetic fields, the differential transmission signal is 
much smaller in comparison to Fig. [sj^a) where spin con- 
serving relaxation taked place. However, due to the ~ B'^ 
enhancement of the spin flipping rate, at larger magnetic 
field spin blockade is suppressed. This in contrast with 
transport experiments where the application of an exter- 



nal magnetic field suppresses the singlet-triplet mixing 
and thus enhances the spin blockade effectP^. 



V. CONCLUSIONS 

We have developed a model describing the trion and 
population dynamics in a photoexcited quantum dot in 
a pump-probe setup. We have included the role of inter- 
subband relaxation including spin flipping and separated 
its role from the spin conserving mechanism. The long 
timescale of intraband spin relaxation leads to a signa- 
ture in the differential transmission signal that is anal- 
ogous to optical spin blockade. In the presence of an 
external magnetic field, the enhancement of spin-flipping 
relaxation rate leads to lifting of spin blockade at shorter 
time scales. This mechanism opens new possibilites for 
the study of spin decoherence processes in semiconductor 
quantum dots with optical probes. 
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Appendix A: Equations of motion 



Here we write the set of equations, derived from Eqs. 



( 12 ) and ( 14) expanded in increasing orders of the optical 



field. Keeping terms up to first order, Eq. (12) for the 
polarization becomes 



)Pi^)-d^(t)(l-<J, (Al) 



where Ena = E-^^^Et 



V^l^ + Uns is the trion energy 



and = 7p + ^np{^cKa + ^sKa + r^)/2 describes the 
trion relaxation rate, which for the p-shell is enhanced 
by the intraband spin-conserving and spin-flipping relax- 
ation terms. 

In second order in the optical field, the equations of 
motion for the electron and hole populations are 
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Finally, for the polarization in third order, we obtain 

ipf; = ^E^^ _ i^_)p(3) + dE{t) [ns^ + ns\ 



5 W'*'^^ 



The last term in the above equation describes contribu- 
tions from interlevel relaxation of electronic populations, 
which as discussed in section IV, leads to spin-dependent 
signatures in the differential transmission signal. 



Appendix B: Analytical Solutions 

The solution for the electronic populations has the 
form 



a' 

X e--^- *0(T)0(t) + e-T- (*+")e(-T)e(t + r)" 
+^(-r)^(t + r)(e-^- (^+") - e-<-'^'+"))] } (Bl) 



where 7^^^ = 7iv + ^npCz^n^rc + Ka^^)- and a^^/ = 
r^jT /(z^sa'Tc + ^scr'Ts)- Comparing the above expression 
with the solution for the hole populations, Eq. (17), there 



are additional terms (ex ^^s) that describe the creation 
of electronic population in the s-shell due to intraband 
relaxation. 

The solution for the third order terms PnJ , which con- 
tribute to the differential transmission signal, is given by 

(for Vnm = 0) 



[\d\\e - 
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for r > and 
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for r < 0. Given that 7^cr ^ Tt?^^? it is clear from the 
above solution that the differential transmission signal 
decays fast for r < 0, while for r > it is dominated by 
the ~ e~^^^^ term at long timescales, leading to spin- 
dependent decay. 
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